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Deterministic many-resonator W-entanglement of nearly arbitrary microwave states 

via attractive Bose-Hubbard simulation 

A. A. Gangat/'0I. P. McCulloch,i and G. J. Milburni 

^ARC Centre for Engineered Quantum Systems, School of Mathematics and Physics, 
The University of Queensland, St Lucia, QLD 4072, Australia 

The Bose-Hubbard model with attractive interactions (ABH) is theoretically known to have a 
phase transition wherein the bosons change from a superfluid phase to a highly entangled nonlocal 
superposition, but observation of this phase transition has remained out of experimental reach. At 
the same time, discrete variable multipartite entanglement of large numbers of resonators is of both 
fundamental and applied interest, but there are no feasible proposals to achieve it. In this theoretical 
work, we jointly address these two problems by (1) proposing an experimentally accessible quantum 
simulation of the ABH phase transition in an array of superconducting circuit microwave resonators 
and (2) incorporating the simulation into a highly scalable protocol that takes any microwave res- 
onator state with negligible occupation of number states |0) and |1) and nonlocally superposes it 
across the whole array of resonators. The large-scale multipartite entanglement produced by the 
protocol is of the W-type, which is well-known for its robustness. The protocol utilizes the ABH 
phase transition to generate the multipartite entanglement of all of the resonators in parallel, and 
is therefore deterministic and permits an increase in resonator number without increase in protocol 
complexity; the number of resonators is limited instead by system characteristics such as resonator 
frequency disorder and inter-resonator coupling strength. Only one local and two global controls 
are required for the protocol. We numerically demonstrate the protocol with realistic system pa- 
rameters, and estimate that current experimental capabilities can realize the protocol with high 
fidelity for greater than 40 resonators. Because superconducting circuit microwave resonators are 
capable of interfacing with other devices and platforms such as qubits, mechanical resonators, and 
(potentially) optical fields, this proposal provides a route toward large-scale W-type entanglement 
in those systems as well. 
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I. INTRODUCTION 

In the standard Bose-Hubbard model [T, repulsive 
on-site interactions compete with inter-site hopping to 
give rise to the well-known superfluid to Mott insu- 
lator quantum phase transition. For integer ratios of 
boson number to total number of lattice sites, the in- 
sulating phase arises when the on-site interaction en- 
ergy is sufficiently larger than the hopping energy so 
that the bosons are localized to single sites and num- 
ber fluctuations on each site are suppressed to zero. The 
superfluid to Mott insulator phase transition was ob- 
served in the landmark experiment of [2 . In contrast, 
the attractive Bose-Hubbard model (ABH) [3Hl] has at- 
tractive on-site interactions and supports a (quasi-) [8] 
quantum phase transition wherein, as the attractive in- 
teractions increasingly dominate over the hopping, the 
ground state |^^s) changes from the superfluid phase 
to a phase in which the bosons are collectively con- 
fined to the same site but nonlocally superposed over all 
sites and number fluctuations at each site are amplified: 
1^..) ^ ;;^(|A^)i|0)2|0)3...|0)M + |0)i|iV)2|0)3...|0)M + 
|0)i|0)2|Ar)3...|0)M + ... + |0)i|0)2|0)3...|iV)M), where Mis 
the number of sites and N is the number of bosons. As 
the ratio of hopping energy to interaction energy goes 
to zero, the ground state becomes exactly a W-state. 



The ABH phase transition has been theoretically consid- 
ered within the context of cold atoms [3l [4] , trapped ions 
[7], and polaritons in cavity arrays [9], but experimental 
demonstration in these platforms remains a difficult and 
unmet challenge. 

Meanwhile, on-chip superconducting circuit systems, 
wherein microwave quanta can be made to interact, have 
emerged as a very effective platform for quantum electro- 
dynamics p!QUT2] . Theoretical and experimental activity 
regarding such circuit-QED systems has recently begun 
to move toward the many-body regime, where quantum 
coherence can potentially be achieved over 1000 or more 
interconnected microwave resonators [12 . This makes 
circuit-QED an attractive option for quantum simulation 
of many-body Hamiltonians. In this theoretical work, 
we show that recent experimental advances can be used 
to realize a circuit-QED system wherein the ABH phase 
transition may be simulated in an array of microwave 
resonators and, unlike in other platforms, the hopping 
energy may be tuned to zero to reach the perfect W- 
state regime. Further, we present a protocol that uses 
the ABH phase transition to convert almost any input 
state It/jin) of a single microwave resonator into a W-type 
state 
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(|V^zn)l|0)2|0)3...|0)M + |0)i|V^,n)2|0)3...|0)M 
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|0)l|0)2|V^in)3...|0)M + ... + |0)i|0)2|0)3...|V^in)M) that 

spans all M resonators of the array. The only funda- 
mental restriction on IV^in) in our protocol is that it have 
negligible occupation of Fock states |0) and |1). The scal- 
ability of our protocol to arrays with large M is due to 
the fact that the ABH phase transition entangles all of 
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FIG. 1. (Color) System schematic. For simplicity, ground connections are not shown. CPW resonators (black shaded cylinders) 
form a one dimensional array. Each resonator has a SQUID (green shaded boxes) connected to its center conductor and is 
coupled to its nearest neighbors with flux tunable couplers (line hatched boxes). The SQUIDS induce a flux tunable nonlinearity 
into the resonator modes. Two bias lines are used to globally tune the SQUIDS and couplers separately. The first resonator 
(not shown) has a separate bias line for its SQUID. 



the resonators simultaneously, rather than one by one, so 
that only one local and two global controls are required 
for the protocol regardless of the number of resonators. 
We estimate that existing technology makes the protocol 
feasible for upto M ^ 40 resonators. Our protocol com- 
plements well the experimental capability demonstrated 
in \TT of producing arbitrary microwave states |V^in) in 
single on-chip microwave resonators. 

The ability to deterministically create many-resonator 
W-entanglement of nearly arbitrary microwave 
states in superconducting circuits is of significance 
for many reasons. For instance, the qubit W- 
state ^(|1)i|0)2|0)3...|0)m + |0)i|1)2|0)3...|0)m + 

|0)i|0)2|1)3...|0)m + ... + |0)i|0)2|0)3...|1)m) is known for 
the robustness of its entanglement under loss [Mj [15] . In 
particular, the global entanglement decay of such a state 
under both phase and amplitude damping are known to 
be independent of the number of qubits [15] . The states 
that are output by our protocol can be understood as su- 
perpositions of qubit W-states, and they therefore have 
the same robustness of entanglement. This is important 
for tests of nonlocality that utilize large scale W-states 
[T6l [TT] . Also, previous proposals for W-state creation 
in superconducting circuits have mainly been restricted 
to qubit degrees of freedom as opposed to the linear 
resonator degrees of freedom that we consider here. An 
exception is the entanglement protocol experimentally 
demonstrated in [l^ for the case of two microwave 
resonators. However, although that protocol may be 
extended to more than two resonators to enable W-state 
creation, it is not feasibly scalable to the many-resonator 
regime. Another significant aspect of our protocol is 
that when li/jin) are coherent states, the protocol outputs 
entangled coherent states (ECSs) [19l[20], which are of 
both fundamental and applied significance in their own 
right [2Qrr24j . For example, coherent states of resonators 
are quasi-classical states, and many-resonator ECSs 
therefore constitute large-scale Schrodinger cats. To our 
knowledge, no feasible proposals for ECS generation in 



superconducting circuits have been put forth prior to 
this, and many-partite ECS generation has not been 
feasibly considered in any platform. Finally, microwave 
resonators in superconducting circuits may interface 
with other types of systems [25^, such as various types 
of circuit and non-circuit qubits [25], mechanical res- 
onators [26, 27 , and (potentially) optical fields [28H35] . 
The protocol presented here may therefore provide an 
indirect route toward large-scale W-entanglement in 
those systems that may be more feasible than other 
more direct approaches. 



II. SYSTEM AND MODEL 

As schematically illustrated in Fig. (IT]), we consider a 
one-dimensional lattice of microwave frequency supercon- 
ducting CPW resonators, each embedded with a SQUID 
that is connected to the center conductor line. The res- 
onators have nearest neighbour coupling via the tunable 
coupler device demonstrated in ^H |37] . As theoretically 
analysed in ^ , the SQUIDS can induce a negative Kerr 
nonlinearity in the microwave modes of each CPW such 
that, looking at only the fundamental mode c of each 
resonator and considering uniform parameters across the 
lattice, the Hamiltonian of the system in the interaction 
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Phases of the ABH as a function of the parameter 
N — ^ ^, c^ Q 25 is approximately independent of 
L function of M as delineated 
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FIG. 3. A plot of Eq. (|3} showing the dependence of the 
ABH parameter r2 on the total number of lattice sites M for 
M> 5. 
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FIG. 4. (Color) Approximate phase diagram corresponding 
to the ABH of Eq. (fl]). The intermediate phase is not shown. 
Dashed blue arrows indicate the last three steps of the entan- 
glement protocol. 



where x is positive ^38j, n may be positive or negative 
[36l [37] , M is the total number of resonators in the lat- 
tice, j is the resonator index, and periodic boundary 
conditions are assumed. Both x ^^d n may be tuned 
in- situ via flux biases, and we designate their ranges as 



< X < X. 



and 



_ ^ _ ^71 



(We note that 



a similar circuit-QED system with open boundary con- 
ditions and no tunable coupling is theoretically studied 
in the driven-dissipative regime in [39^.) Experimental 
considerations related to damping, higher order terms 
of the nonlinearity, and resonator frequency disorder are 
discussed in Section W\ 

Equation (fl]) is precisely the attractive Bose-Hubbard 
model (ABH) that is theoretically studied in [3]-[7]. For 
a fixed total number of quanta N and A^ > 1, the ABH 
phase changes qualitatively as a function of the param- 



eter r = 
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as depicted in Fig. (2). ri 
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X(A^-I)' 

is a characteristic constant of the ABH that is largely 
independent of lattice size [4l|6]). When r < ri the at- 
tractive on-site energy dominates the hopping energy and 
the quanta in the eigenstates are collectively localised to 
single sites but superposed across all sites to form the 
nearly-degenerate W-states: 
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,i/c(27rj7M) 



\N)j n |0)- (2) 
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where k is an integer in the range < /c < M — 1. These 
W-states become the exact eigenstates as r ^ 0, and our 
proposal for ABH simulation is the only one to date that 
is able to access this perfect W-state regime as it is the 
only one in which k. is tunable to zero so that r = 0. r2 is 
another characteristic constant of the ABH such that for 
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FIG. 5. (Color) Approximate phase diagram corresponding 
to the single-site nonlinearity model of Eq. Q. Dashed blue 
arrows indicate the first three steps of the entanglement pro- 
tocol. 



T > T2 the ground state is a superfluid state. It is always 
the case that ri <T2- When n < r < r2 the ABH phase 
is intermediate in nature between the W-state phase and 
the superfluid phase. T2 increases monotonically with M: 
T2 = Ti for M = 2 [6 , Ti<T2< 0.3 for 3 < M < 5 [5], 
and 

T2 ^ [2Msin^(iT/M)]-^ (3) 

for M > 5 [4 , which is approximately linear in M and 



is plotted in Fig. (|3|. The approximate phase diagram 
for the ABH is represented in Fig. Q with the regime 
Ti > r > r2 omitted. 

We also consider the case when only the first site in the 
lattice may have a non-zero attractive interaction with 
strength xi and range < Xi < Xmax- 



number of quanta N in the system. \il)in) therefore has 
the constraint that |Cn| ^ for n < 2. As mentioned at 
the end of section [TTj each total quanta number n from 
the distribution of N may be treated independently. We 
may therefore express the total system state at all times 
as 



H 



M 



For A/" > 1, the two parameter regimes of interest are 
Xi{N — l)/2 ^ \i<i\^ in which the quanta in the lowest 
energy eigenstate are localized only at the first site, and 
Xi(A^ — l)/2 <C |a>:|, in which the quanta in the lowest 
energy eigenstate form a superluid state across the whole 
lattice for any finite value of n. The approximate phase 
diagram for this model is represented in Fig. ([5|. 

Through local control of the nonlinearity Xi of the first 
resonator and global control of Xj^ and n^ our proposed 
system can simulate both Eqs. (IT]) and Q. An initial 
state of the proposed system with A" > 2 in the low- 
est energy eigenstate and zero quanta in the other eigen- 
states may therefore adiabatically transition between the 
single-site localized phase, the delocalized phase, and the 
W-state. If instead the system starts in a superposition 
of different A^ in the lowest energy eigenstate, the linear- 
ity of quantum mechanics dictates that each component 
of the superposition with N > 2 may independently un- 
dergo the phase transitions. The entanglement protocol 
revolves around these phase transitions, and therefore 
necessarily excludes initial states with number distribu- 
tions that have significant occupation of A" < 2. 



III. ENTANGLEMENT PROTOCOL 

Mention is made in the literature of using adiabatic 
transitions of Hamiltonians to generate entangled Fock 
states [9] HOI [41]. In particular, Hartmann et al. [9 sug- 
gest using an atom-cavity realization of the ABH to cre- 
ate a polaritonic (photon-atom hybrid) W-state of Fock 
states in multiple optical cavities via adiabatic transition 
from the superfluid regime (r > r2) to the W-state regime 
(r < Ti) by tuning x- Here we show that the ABH can be 
used to create W-type entanglement of nearly arbitrary 
pure photon states. 

With the system in the vacuum state and x = ^ = 0^ 
the first resonator of the system is initialized in a state 
IV^in) = Cn\n)^ where \n) denotes a Fock state with n 
quanta and Cn are complex amplitudes, so that the initial 
total system state is 



M 



l*(0)) = |V'm)i=in®|0)i- 
i=2 



(5) 



The probability distribution of |V^in) in the number ba- 
sis is therefore the probability distribution of the total 



i*w) = E^»i*"W)' 



(6) 



n=0 



where |^n(t)) denotes the system state at time t with def- 
inite total quanta number n, and the evolution of each 
|^n(^)) niay be considered separately. The entanglement 
protocol employs manipulations of Xi, X^ and hz that si- 
multaneously evolve each |^^(t)) from |^^(0)) = \n)j=i 
(a Fock state of the first resonator), into |^n(^)) = 



^M 



/M 



XI 7=1 \^)j Ylry^j |0)r (a W-state spanning all of the 
resonators), for all n > 1. The total system state at the 
end of the protocol is therefore a W-type distribution of 
the input state Itpin)- 
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(7) 
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Because the entanglement protocol maintains each 
|^n(^)) in a single eigenstate of the system, we may ex- 
press the time-evolution as 



l*n(i)) 



-i<Pnit) 



n{t)) 



eig7 



(8) 



where ^n(^) = \ Jq ^n{t')dt' and En(t) denotes the 
energy of the eigenstate |n(t))eig that is occupied by 
the n quanta at time t. We first describe the tra- 
jectory of |n(t))eig during the entanglement protocol 
to explain how it evolves from |n(0))eig = |^)j=i to 
|n(T))eig = -^ Ej^i W)j Yir^j |0)r, and then explain 
how e-^^-^^) = 1 is achieved so that \^n(T)) = 
-^ T.f=i \^)j Ylr^j \^)r- In general (t)n{t) 7^ (/)^(t) when 
n y^ m^ which can result in a distortion of the phase 
information contained in the initial state. 

The steps of the entanglement protocol are illustrated 
with numbered arrows in the phase diagrams of Figs. J4]) 
and ([5| and in the pictorial representation of Fig. (I6|i, 
which depicts the evolution of the full system state |^(t)) 
(see Eq. (Ig])). Let tx denote the time in the entanglement 
protocol after step number x. The first three steps of the 
protocol (see Figs. ([5| and (|6|) use the model of Eq. 
Q to convert each |n) localized in the first resonator 
into a delocalized superfluid state across the whole lat- 
tice. In the first step, Xi is rapidly increased to Xmax in a 
time Ati, which does not change the occupied eigenstate: 
|^(^i))eig = |^(0))eig- The attractive interaction thus in- 
troduced into the first resonator provides an energy bar- 
rier against which k may be tuned adiabatically in step 2 
from to -f^max inatime At2 > 4:K.max/{Xmax{n-l))'^. 
In step 3 xi is adiabatically tuned to zero in a time 
At3 > Xmax{N - \)I2kI^^^. Duc to the negative sign 
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FIG. 6. (Color) Illustration of the entanglement protocol for the case of four resonators. Square boxes denote single resonators, 
full black lines shared between square boxes denote inter-resonator coupling [n) off, dashed black lines denote coupling on, 
green shading denotes resonator nonlinearity (x) on, white shading denotes nonlinearity off, periodic boundary conditions are 
implied, and dashed blue arrows denote steps of the entanglement protocol. \il^ln) — e~^'^'^'^^^Cn\n) denotes a version of \il^in) 
with modifications to the phases of each number component. After the first step, \il^in) may adiabatically transition to the 
superfiuid state by turning n on (step 2) then turning xi off (step 3). After rapidly switching the sign of k in step 4, IV^J^) is in 
the superfiuid groundstate. Step 5 (x is adiabatically turned on across the whole lattice) and step 6 [n is adiabatically turned 
off) induce the ABH phase transition so that \^|J*n) becomes nonlocally superposed in each resonator (W-type entanglement). 
Finally, turning off the resonator nonlinearities x (step 7) restores the lattice to an array of uncoupled linear resonators that 
now contains the entangled state. With correct timing of the last step, (j)n{T) = (T is the total time of the protocol) and 
IV^*^) is restored to \'ij^in)- 



place |n(t5))eig in any of the three phases of the ABH (see 
Fig. ([2|) depending on the exact values of Xmaxi i^raax-, 
and n. In step 6 a^: is adiabatically tuned from Umax to 
zero in a time Ate (see section |VB for the adiabaticity 



constraint) so that r = and |n(t6))eig is the W-state 



of K. and the single-site nonlinearity during the first three 
steps, |n(t3))eig is actually the highest energy eigenstate 
of the system. In step 4 a^: is rapidly tuned from —Umax 
to -^nmax in a time At4 so that |n(t4))eig is the ground- 
state of the system. I n(t4))eig can therefore adiabatically 

connect to the lowest energy W-state |^y^~ ) of Eq. (l2|. \^w~ )• ^^ ^^^^ step, although the eigenstates become 
Steps 5 and 6 of the protocol (see Figs. Q and (p|) ac- 
complish this by using the ABH (Eq. ([TJJT^ In step 5 the 
global parameter x is adiabatically tuned from zero to 
in a time Ats (see section 



X' 



ity constraint) so that r 



VB| for the adiabatic- 
r:/Xmax{n-l). TMs may 



degenerate as r ^ 0, the on-site attractive interaction 
energy proportional to x serves as an energy barrier that 
protects the system from transitioning out of the ground 
state as long as the tuning is adiabatic with respect to 
X- As /^ is now off, the superposition is locked into place 



and X ^^y be rapidly tuned in step 7 from Xmax to zero 
in a time Atr without altering the occupied eigenstate: 
|n(T))eig = -^ Y.f^i W)j Yir^j |0)r. As the nonlinear- 
ity is now off, the CPWs are not hybridised with the 
SQUIDS and the state of the system is purely photonic. 

The effects of the entanglement protocol on (pniT) may 
be understood by separately considering the contribu- 
tions of the on-site nonlinearity terms and the nearest 
neighbour hopping terms of Eqs. ^ and Q. In the case 
of a single site, the nonlinearity contributes a number 
dependent phase n J^ ^^^ ^ dt' J42] so that appropriate 

timing can make /^ ^^"^ ^ dt' equal to an integer multiple 
of 27r and the effective contribution to 0n(^) is zero for 
all n. It may be conjectured that, with appropriate tim- 
ing, the same type of cancellation for all n could occur in 
the multi-site case. However, this is a far less trivial con- 
sideration due to the changing distribution of the quanta 
with time. Remarkably, numerical investigations show 
the conjecture to be exactly correct. For the hopping 
terms, it can also be conjectured that the contribution to 
(j)n{T) can be approximately cancelled by simply making 
At2 = Ate and Ats = Ats so that the phase evolution 
due to the hopping terms when n <{) (steps 2 and 3) may 
cancel the phase evolution due to the hopping terms when 
/^ > (steps 5 and 6) because of the opposite signs of the 
hopping energy and because step 5 is in some sense the 
reverse of step 3 superposed on each site. Remarkably, 
numerical investigations reveal that this cancellation is 
also exact, even in the absence of adiabaticity. It is there- 
fore possible to achieve (pniT) = through appropriate 
timing of the steps of the entanglement protocol, as we 
demonstrate in the numerical simulations below. As an 
alternative to the cancellation of (j)n{T) through timing, 
a calibration may be done whereby the phase (pniT) is 
measured after test runs of the protocol with {ipin) = |^)i 
for different n so that the arbitrary input states gener- 
ated by, for example, the method in [13 may be prepared 
with appropriate offset phases for each number compo- 
nent so that the (j)n{T) accumulated through the protocol 
are all cancelled. 



EXPERIMENTAL CONSTRAINTS 



IV. DETECTION 



Verification of entanglement may be done by employ- 
ing bipartite Wigner tomography [18 between different 
pairs of resonators to reconstruct their joint density ma- 
trix. Alternatively, if after the entangled state is prepared 
the system is made into a linear network by quenching 
K. from zero to ^maxi the protocol of Tuferelli et al. [43] 
allows for a single qubit tunably coupled to any resonator 
to be used to reconstruct the state of the entire array. 



With present-day devices [36H38] , Kmaxl'^^ and 
Xmax I '^^ may reach as high as hundreds of MHz and 
may be tuned on a timescale of a few nanoseconds. How- 
ever, K should be limited to about nmax/"^^ = 30 MHz 
in the interest of preserving high experimental fidelity of 
the hopping Hamiltonian with the tunable coupler [44^ . 
Modes Cj have fiux-dependent Ti and T^ due to fiux- 
dependent hybridization of the CPWs with the SQUIDs. 
Q factors of over one million have been demonstrated for 
on-chip CPWs [45 . This gives Ti > 20 /is for the 27r x 7.5 
GHz CPWs that we assume here. As mentioned in J38] . 
the SQUIDS embedded in the CPW resonators should be 
able to achieve very long coherence times considering the 
recent Ti and T2 measurements of the Josephson junction 
qubits in [46ti48] . The 2D Xmon qubit demonstrated in 
[48 , for example, shows a maximum Ti of approximately 
44 /is and Ti ^ T0/2 ^ 20 /is at the fiux insensitive point. 
Further, considering an asymmetric SQUID as in [38 al- 
lows for T(j) > 1 ms at the fiux value where x is maximum 
[49 . For our purposes here we therefore approximate a 
constant Ti value of 20 /is and a fiux-dependent T^ value 
that varies linearly from 1 s to 300 /is as x is increased 
linearly from to Xmax- 

In this section we show that these and other experi- 
mental capabilities make our protocol feasible for array 
sizes of greater than 40 resonators. In particular, for 
the choices of Kmaxl'^^ = 30 MHz and Xmaxf^^r = 14 
MHz, the ranges 2 < n < 40 and 2 < M < 42 become 
available, which can accommodate high amplitude many- 
partite ECSs. 



A. Higher order terms of CPW nonlinearity 

The full expression for the nonlinearity {Hnl) in- 
troduced into the CPWs by the embedded SQUIDS is 
given in Eq. (27) of [38 . Looking only at the fun- 
damental mode, and using the relations <l>o = /i/2e, 
Ej = il^f/Lj, Lj = L'Inu and L' = 2E'c / {ul^e') 
we find 



^^^ = |:^(^(^:;^j rii^cAc'^c] , (9) 

where ^z (0 ^ t^z ^ 1) is the inductive participation ratio 
dependent upon the fiux through the SQUID loop and 
uJc,m is the frequency of CPW mode m in the absence of 
the SQUID. After the rotating wave approximation, to 
be able to neglect nonlinear terms higher than {c^cY we 
require 



(^X^Ag/ A' » l(^^)'20(n_.)', (10) 
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(11) 



Larger values of Xmax come at a cost of reduced n^nax 
but allow for quicker adiabatic tuning of n as per the 
discussion below. Smaller values of Xmax enable larger 
array sizes, also as per below. Larger values of Umax 
may be accommodated for a fixed Xmax by using higher 
order CPW resonator modes. If we assume a funda- 
mental mode frequency of co'c,i/27r = 7.5 GHz, selecting 
Xmax I '^^ — 14 MHz allows Umax = 40 whcu the funda- 
mental mode is used, selecting Xmax I "^^ = 25 MHz allows 
^max = 23 when the fundamental mode is used, and se- 
lecting Xmax I '^^ — 100 MHz allows Tiraax = 6 whcu the 
fundamental mode is used. 



B. Adiabaticity 

Fig. (7) of [4 reveals a universal behavior of the ABH 
in which the population of the fundamental normal mode 
of the lattice decreases approximately linearly with the 
parameter r/r2 in the region < r/r2 < 1, and is ap- 
proximately constant in the region r/r2 > 1. This trans- 
lates to the two adiabatic constraints ^^^ <C Xmax and 

-^-^ <C 2Kmax in the respective regions. 
We denote the value of r at the end of step 5 as 
-. The consideration in the next subsec- 



Xmax(n-l) • 

tion shows that it is of interest to minimise the time taken 
to tune r from T2 to t\. We find that this occurs when 
r* ^ T2 (although this is not simultaneously possible for 
all n, it is sufficient to assume so as a rough estimate). 
In this case the adiabatic constraint to tune x from 
to Xmax in a time Ats is ^ < \x^^2{n - 1). Since the 
constraint only becomes applicable as r approaches r2, 
the approximation x ~^ x" is made on the RHS to yield 



A^5 ^ ^/ f^max- 

For step 6 we have the requirement ^ 



(12) 



< — y^ (n 

- l/imaxy'^ 



l)r2. Tuning A>: from /^„ 



dt ^ 10 ^max\ 

to i^{ri) in a time At^a requires 



AtQa ^ 10 ^^ ^^ , and tuning /^ from K:{ri) to zero in a 
time Ateb requires Ate^ > 10 — — — . The constraint on 



step 6 is therefore 



Ate > 10/xr 



(13) 



For Xmax/'^TT and Kmaxl"^^ on the scale of tens of 
MHz, this puts steps 5 and 6 on the 10-100 nanosecond 
timescale, which is well below the estimated T\ and T^ 
of the resonators. 



C. Disorder 

Disorder in inter-resonator coupling frequencies can be 
neglected when resonator frequencies greatly exceed the 
coupling frequencies [50 , which is the case here. We 
focus therefore upon disorder in resonator frequencies. 

The inevitable frequency spread Auj of the resonators 
relates to the discussion in [3 concerning the fact 



that non-degeneracy of the lattice sites induces non- 
uniformity in |eig(n,t)) in the site basis at a rate 1/Acj. 
This is relevant in the intermediate regime t\ < t < T2^ 
which is traversed in step 6 given the assumptions above. 
If this intermediate regime is traversed on a time scale 
comparable to or longer than l/Acj, the quanta have 
time to localise in the lowest energy site. In order to 
achieve symmetric W-entanglement with high fidelity, 
we therefore require Atint ^ l/Acj, where Atint is 
the time it takes to traverse the intermediate regime. 
From the previous subsection, Atint = At^a- Letting 
AUnt = l/lOAo;, we find -£ > ^{1 - ri/r2), where 
it is assumed that r2 = ^maa. ^ ^j^^g g^ygg 



Auj< 



vmax'^max 



im{Kr. 



Xmax y^' 



max y^rmn 



l)/4)' 



(14) 



which gives a lower bound on Xmax- Smaller values of 
Xmax (oi" ^min) enable larger r2 (and therefore larger M) 
for a fixed rimin (or Xmax)^ but also require smaller values 
of Auj and larger At^. Current capabilities demonstrate 
Auj > 27r X 1 MHz for GHz frequency CPW resonators 
[bOy Choosing rimin = 2 and hZraaxf^^r = 30 MHz, we 
thereby find the constraint Xmaxf^^r > 14 MHz. Using 
Eqs. ([3| and (11) this translates to M < 42 and rimax ^ 
40, which is sufficient for many-partite ECS creation. If 
Auj is reduced by a factor of two, we find Xmox/27r > 
7.5 MHz and M < 80 for rimin = 2. Reducing Auj to 
27r X 0.1 MHz gives Xmax/'^Tr > 1.6 MHz and M < 360 

lOr Tijjiin — ^' 



VI. SIMULATIONS 

We demonstrate the entanglement protocol by numer- 
ically integrating the master equation 



dp 
Jt 



M M 

[H, p]+J2 7fr'D[cj]p + J2 7^S[cj]p, 



(15) 



where p is the density matrix for the CPW chain, H is the 
system Hamiltonian, V[cj]p = CjpCj — CjCjp/2 — pCjCj/2 
is the amplitude damping operator, G[cj]p = c^-Cjpc^-Cj — 
{cjCj)'^p/2 — p{cjCj)'^/2 is the phase damping operator, 
1/Ti is the amplitude damping rate, and l/T^ is the 
phase damping rate. The integration is done using the 
fourth order Runge-Kutta method with a timestep of size 
10~^^ s. As discussed in section|v| we model the damping 
with a constant Ti value of 20 ps and a fiux-dependent 
T^ value that varies linearly from 1 s to 300 ps as x is 
increased linearly from to Xmax- 

In the first simulation, a system of three resonators 
(M = 3) is considered with an input state IV^in) = 
^|2) + 7^ 1 3). The system parameters Xmax I "^^ = 100 
MHz and Kmaxl'^^ = 30 MHz are used, which satis- 
fies the constraints discussed in section [V| and the tim- 
ings of the protocol steps are chosen to respect the adi- 
abaticity constraints. The Hilbert space at each site 
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FIG. 7. Fidelity of the system state |^(t)) with the tar- 
get state l^^w) = -^Y^^j=i\'4^in)jY[r^j^\0)r during the fi- 
nal step of the entanglement protocol with system param- 
eters Xmax/27V = 100 MHz, Kmax/^TV = 30 MHz, M = 3, 

and damping as explained in the text. The input state is 
iV^in) = ^1^) "1" ^1^)- ^^^ fidelity oscillations are due to the 
phases 02 (t) and 03 (t) of the respective number components 
evolving at different rates due to the number-dependent fre- 
quency induced by the Kerr nonlinear ity at each site. The 
frequency of the oscillations decays to zero as x is tuned to 
zero, and the timing Atr of the step is chosen such that the 
oscillations cease at a fidelity peak of about 97.5%. 



is truncated at n = 3. The fidelity |(^(t)|^vF)P of 
the system state |^(t)) with the target state l^^w) = 
75 Yfj=i \'^in)j Ylry^j ^|0)r during the final step (step 7) 
of the protocol is shown in Fig. ^. The total physical 
time of the simulation is 0.1064 jas. The same simula- 
tion is also performed without damping (not shown) and 
a peak fidelity of 98.6% is found. Further numerical in- 
vestigations (not shown) reveal that, in the case of no 
damping, the peak fidelity is limited only by imperfect 
adiabaticity. 

In the second simulation we use M = 3, Xmax/27r = 
120 MHz, and K.max/^7r = 30 MHz. The input state 
is 1^,^) = ;7g|2) + ^|3) + e'"/^^|4), and the Hilbert 
space at each site is truncated at n = 4. The fidelity 
|(^(t)|^iy)P during step 7 is shown in Fig. ^. The 
timings of the protocol steps are based on trial simula- 
tions that determine the level of adiabaticity needed to 
achieve the high peak fidelity shown. The total physical 
time of the second simulation is 0.3556 /is. 



VII. CONCLUSIONS AND OUTLOOK 



FIG. 8. Fidelity of the system state |^(t)) with the tar- 
get state I'^w) = -^Y^^j^i \'^in)jYlr^j^\0)r during the fi- 
nal step of the entanglement protocol with system param- 
eters Xmax/27T = 120 MHz, Kmax/^TT = 30 MHz, M = 3, 

and damping as explained in the text. The input state is 



I^^^) = 7ll2) + ^ 



|3) + e^"/^^|4). The fidelity oscillations 



are due to the phases of the number components evolving 
at different rates due to the number- dependent frequency in- 
duced by the Kerr nonlinearity at each site. Each period of 
the oscillations has three peaks due to the three frequency 
differences of the three number components of iV^in)- The 
frequency of the oscillations decays to zero as x is tuned to 
zero, and the timing Atr of the step is chosen such that the 
oscillations cease at a fidelity peak of about 95%. 



in-situ tuning of the hopping energy to enable access to 
the perfect W-state regime. We have further presented a 
protocol built around the attractive Bose-Hubbard simu- 
lation that deterministically produces W-type entangle- 
ment of nearly arbitrary single resonator states over a 
very large number of microwave resonators. We have nu- 
merically demonstrated our protocol with complex input 
states in an array of three resonators using realistic pa- 
rameters and shown the attainability of high output state 
fidelity. The highly entangled large-scale states produced 
by the protocol have both fundamental and applied sig- 
nificance. 

Looking ahead, this protocol may also offer a path to- 
ward creating large scale W-type entanglement in qubit, 
mechanical, and optical degrees of freedom due to the po- 
tential that superconducting circuits hold for interfacing 
with those platforms. It is also interesting to consider 
if improved scaling of resonator number with system pa- 
rameters could be achieved for lattice geometries other 
than the one-dimensional case considered here. 



We have proposed a cicuit-QED simulation of the at- 
tractive Bose-Hubbard model in which the as-yet ex- 
perimentally unobserved superfluid to W-state quantum 
phase transition may be realised with existing technol- 
ogy. A unique aspect of this proposal is the capability of 
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